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00 | Abstract 

We consider parabolic partial differential equations of Lotka-Volterra type, with a non-local 
nonlinear term. This models, at the population level, the darwinian evolution of a population; 
Oh I the Laplace term represents mutations and the nonlinear birth/death term represents competition 

■ leading to selection. 

Once rescaled with a small diffusion, we prove that the solutions converge to a moving Dirac 
mass. The velocity and weights cannot be obtained by a simple expression, e.g., an ordinary 
differential equation. We show that they are given by a constrained Hamilton- Jacobi equation. 
This extends several earlier results to the parabolic case and to general nonlinearities. Technical 
new ingredients are a BV estimate in time on the non-local nonlinearity, a characterization of 
the concentration point (in a monomorphic situation) and, surprisingly, some counter-examples 
showing that jumps on the Dirac locations are indeed possible. 
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1 Introduction 



- T— I ■ 

This paper is devoted to study the asymptotic behavior of solutions to parabolic Lotka-Volterra 
■ equations. They describe the dynamics of a population density n(t, x) which expands (or decays) 



with a growth rate R(x,I{t)) which changes sign. In the theory of adaptive evolution ([27, 
[26j [17] lllj). the parameter x denotes a physiological "trait" and I(t) an environmental unknown 
shared by the total population and which is used as a nutrient. The population can use the nutrient 
I(t) differently, depending upon the trait x, which makes that certain fitter traits should emerge, 
modifying the environment and thus allowing other traits to emerge. This mechanism uses mutations 
in the population that we model here by a mere diffusion; but more realistic integral kernels can be 
handled as well, \18\ I12j. This elementary modeling leads to the equations 

|n £ -eAn £ = — R(x,I s (t)), x 6 t>0, 

dt e (1) 

n £ (t = 0) = rt° G L 1 (R d ), n° £ > 0, 
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I e (t) = / ijj(x)n £ (t,x)dx. 



(2) 



The function ip is given and measures the 'predation' of individuals with trait x on the environment. 
Notice that we have rescaled the problem in order to include the idea that mutations are small (or 
rare) thanks to the small parameter e > 0. Such models, together with related asymptotic, can be 
derived form individual based stochastic processes in the limit of large populations; we refer to 1 1 3 [, II 4 j . 
Also the problem without diffusion (a coupled infinite system of elementary differential equations) is 
interesting from the point of view of large time behavior; we expect that the dynamics concentrates 
on large times and several related results can be found in the literature, see |16|. [28J for instance. 

Our purpose is to show that, under various assumptions on the rate R(x,I), the population n e 
concentrates as a Dirac mass (or a sum of Dirac masses), a mathematical way to express that well 
identified species emerge from the adaptive landscape defined by the rate R. The interesting feature 
being to describe the dynamics of these Dirac masses. Here, we perform this analysis in a completely 
rigorous manner. But the idea to analyze adaptive dynamics in those terms, and the formalism, goes 
back to [18j . where a well founded biological system (the chemostat) was studied. In [6], we performed 
a rigorous asymptotic analysis in the case of integral operators and with linear dependency on the 
environmental unknown I(t). The case of a system, for a population with adults and juveniles, was 
studied in [12]. 

Although related, the situation of reaction-diffusion systems, as they arise in combustion, is quite 
different. The simple Fisher-KPP equation may serve as a model, which amounts to use R = 1 — n e 
in (pQ), and is known to lead to the propagation of a front. This means that in the limit e — > 0, the 
solution converges typically to either (uncolonized region) or to 1 (fully colonized region), and the 
transition occurs on the "front", i.e., a hypersurface which dynamics can be described by the level set 
of a solution of a Hamilton- Jacobi Equation (see [5] [7] [8] [20], [29J |22j ) . Therefore, the limiting objects 
are, geometrically at least, very different in the cases of reaction-diffusion equations and of Lotka- 
Volterra equations. This is the reason why in the later case a new type of equation, the constrained 
Hamilton-, Jacobi equation, occurs to focus on the isolated points of the Dirac locations. The theory 
of viscosity solutions to Hamilton-Jacobi equations occurs naturally in our derivation because a phase 
occurs naturally; it is quite elaborate now and general introductory references are [U [15] [3] [19] [21]. 
But a large part of the present paper uses other general ideas that can be read without knowledge of 
this notion. 

The paper is organized as follows. We first state (section [2]) simple and general results and in 
particular the convergence result to a Dirac mass for dimension one in "monomorphic situations". 
These results are completed by BV bounds on I e (t) (section [3]) and by the asymptotic analysis through 
constrained Hamilton-Jacobi equations (section HJ, in any dimension. Consequences are drawn in 
section [5] they imply the concentration of n £ as a sum of Dirac masses in a general multidimensional 
setting. Section [6] is devoted to counterexamples to the continuity of I(t) (the limit of I £ ) and of the 
dirac locations, and to "smallness" conditions implying continuity. Several further results are presented 
afterwards: the monomorphic situation is completed in section [7J the case of several environmental 
variables is treated in section [8] 

2 Assumptions and main results 

In this section, we give assumptions for the coefficients and data arising in ([I])-© and state the main 
convergence result in a particularly simple case. More general results follow from the analysis and 
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proof we perform later; these results, which are more technical, are stated along the paper. 
We assume that there are two constants ip m , ipM such that 

< Vm < V> < i>M < oo, ip G W 2 ' 00 ^). (3) 

The quantity R is called the invasion exponent because it describes the ability of the individuals of 
trait x to invade the population with environmental state lit)- We assume that there are two constants 
< I m < Im < oo such that 

min R(x,I m ) = 0, max R(x, Im) = 0, (4) 
and there exists a constant K > such that, for any i £ R 1 ', I £ I, 

-K<—(x,I)<-K' 1 <0, sup WR^rjW^ai^KK. (5) 

ai I m /2<I<2I M 

We will also use the assumption 

n°GL°°(M d ), Vn° £ eL l (M. d ) and J m < f ip(x) n° E (x)dx < I M , (6) 

Jw. d 

and the notation 

g e (t) = / n £ (t,x)dx. (7) 



Notice that the assumption ([3]) and the bound I m < J Rd tp(x) n £ (t,x) dx < Im imply 

-j— ='■ Pm < / n £ (t,x)dx < p M ■= ~r~- 
WM jRd 1pm 

One can have in mind the particular but more intuitive example 

ij> = 1, R(x, I) = b(x)Q 1 (I) - d(x)Q 2 (I), (8) 
with Qi G C^R) for i = 1, 2, and 

Qi(0<o, Q' 2 (-)>o, g l >o. (9) 

b>b m >0, d>d m >0 and b, d £ W 2 >°°{R d ). (10) 

We recall that we have the following existence result and a priori bounds (here C denotes various 
constants which maybe different from line to line) 

Theorem 2.1 With the assumptions (Ep-|7|) and I m — Ce 2 < I £ (0) < Im + Ce 2 , there is a unique 
solution n £ G C ; L l (U. d )) , to equation (QP-(0), and it satisfies, 

Im-Ce 2 <I £ (t)<I M + Ce 2 . (11) 
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A proof of existence can be found, for instance, in [16] and we do not recall it. The uniform bound 
can also be also be found in (6J [28] (according to [16] the lower bound is unessential because they 
are indirect ways to guarantee non-extinction a posteriori, we keep it here for the sake of simplicity). 
Section [3] gives another and stronger uniform bound in time, a uniform BV bound on I e {t). Here, we 
just indicate the derivation of the upper bound in (|lll) . We have 

I ip(x) n s (t,x) dx = ej Rd n £ (t,x) Aip dx + - / n £ (t, x)R(x, I e (t)) 
at jRd e J R d 

< Cl eI e {t) + -Ie(t) max R ( x >*f» 

and, from the assumption (0), the right hand side becomes negative as soon as I £ (t) overpasses 

Im H £ 2 and the result follows. 

K 

We can state a very simple version of our results in the simple case when dimension d is equal to 1 
and when, typically, the function R(x, I) is monotone in x. 

Theorem 2.2 (Dimension d = 1) We assume (Ep-fSJ), the technical condition on n £ in Theorem 



4-1 below, and 

VT m < I < Im there is a unique X(I) el such that R(X(I),l) = 0. (12) 

Then, the solution n e (t) to equation (Cp-(0) converges in the weak sense of measures (see also Remark 
[X^l below) 

n ^(*) 1^ Q(t)5(x-x(t)), (13) 

and we have 

x(t) = X(I(t)), R{x(t),I(t)) =0, 
and the pair (x(t), I(t)) satisfies the constrained Hamilton- J acobi equation given later on (see Section 

When R(x, I) has the special form insuring uniqueness in Theorem \4-l\ then the full family n £ 
converges. 

Such a population is called monomorphic because a single trait is represented asymptotically. This is 
the general situation with a single environmental variable I(t) (this is called the Competitive Exclusion 
Principle, [32]). To go further, several environmental variables can be introduced (see [181 126] and our 
results in Section [8]). 

3 BV estimates on I £ (t) 

As a first step in our analysis of the limit e — ► in (pQ), we prove strong convergence of I £ {t). Therefore, 
we complete result of Theorem 12.11 by the: 

Theorem 3.1 With the assumptions {3p-(E|), we have additionally to the uniform L 1 bound HI]) , the 
local uniform BV and sub-Lipschitz bounds 



j t h{t) > -eC + e- K ^ J ^ X )n° £ (x)^^- 
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d qS) > -Ct + f (l + iPixMix) 11 ^' 1 ^ 



dt 



where C and are positive constants. 

Consequently the extracted limits satisfy that I(t) is nondecreasing as soon as there exist a constant 
C independent of e such that 



V^K°(x) ^ X,40) > -Ce°«/ £ . (14) 



Then, we also have, for all T > ; 

f [ n £ (t,x)R(x,I e (t)) 2 dx dt <C e (1 + eT). (15) 

Jo JR d 

Remark 3.2 Our condition {1$ on the initial data leaves place for a possible initial layer. In order 
to avoid it, one may choose well-prepared initial data that satisfy R (x, Ig) = (for instance by tuning 
the total mass ofn®). When it is not fulfilled one might observe, e.g. in numerical simulations, a fast 
variation of I e (t) for t rj 0. 

Remark 3.3 From the bound $15\). we can deduce that the weak limit n in L°°(M. + ; M 1 (M d )) of n £k 
is supported (a.e. in t) only at points x such that R[x, J(t)) ■ We give a more precise result below. 

Remark 3.4 Notice also that one can prove following the same lines that the family J Kd n e (t, x)^f(x)dx 
is bounded in BV\ oc for any W 2, °° test-function \I/ . The extracted limits n £k therefore converge a.e. in 
time and not merely in w~L°° (M. + ;w-M 1 ( 



Remark 3.5 The interested reader can check that the proof below also extends previous results of 
JUl f or continuous differential system 



d f 

—n = nR[x,I(t)), I(t) = / n(t,x)dx. 



When R is monotonic in I and assuming there is a unique x such that = R(x, Im) = max a . R(x, Im), 
then as t — > oo, we have n(t) — 1 pm5(x — x) with Im = Pm4>{x). 



Proof. We begin with the proof for I £ , then we show (|15p and finally indicate the variants for the 
proving the result on g e . 

Concerning I £ , we multiply the equation by ip and integrate over R , this yields 

—IS) = £ J n £ (t, x)AiP(x) + JS), (16) 

with J £ {t) defined by 

JS) = J y tP(x) R(x,iS))- 

The integration by parts is justified because ip G W 2 '°° and, for e fixed, both n £ and Vn £ belong to 
L 1 by easy a priori manipulations. 
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The result relies on an estimate on J £ {t). In the same way, we have 

rj 77 9 

-J £ (t) = J n £ A^(x)R(x,I £ (t))]dx + J ^(x)R(x,I E (t)) 2 dx 

(17) 

+ J^(x)g- I R(x,I £ (t))dx -I £ (t). 

Now we use (|16p to recover J £ {t) from —I £ {t) in the last term. And we notice the following properties : 
by ([5]) and (fTT|) . we have 

J n £ A[iP(x)R(x,I £ (t))] - J n £ (t,x)AiP(x) J n £ (t,x)^(x)^R{x, I £ (t)) = O(l) > -K x , 

Therefore, using again ©, 

(1 IT) 9 / 77 (7 

-J £ (t) = 0{l) + J ^(x)R(x,I £ (t)) 2 dx + J -j^{x)—R(xJ £ {t))dxJ £ {t). (18) 
But we have 

J n £ (t,x)iP(x)-?jR(x,I £ (t)) < —K 2 < . 
Using also that the second term in the right-hand side of the above equality is positive, we obtain 

From this differential inequality we find 

(Je(t))_<e^- + (j £ (0))_e- K ^. (19) 
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The local i?V bound on I £ (t) = J ip(x)n £ (t,x)dx, as well as its monotonicity in the limit, follows by 
inserting this inequality in (|16j) . 



In order to prove (I15jl . we ague as follows. We have from (|16p 

/ J £ (t) = I £ (T)-I £ (0)+eO(T) = C + eCT. (20) 
j o 

On the other hand, from (|18p. we also have, for some K(t) > bounded from above and from below 
away from 

±J e (t) = C+ [ ^( x )R{xJ £ {t)fdx J £ (t). 
dt J Rd e 2 e 

After integration, this gives 

J £ (t) = J £ (0)e~ !o ~ + I e ~ Is ~ ( C + J Rd !h ^-^)R(x,Ie(s)) 2 dx^j ds, 



and thus, using assumption (|14p and with K = min K(t) > 0, we arrive at 

[ T Je(t)dt > -C + F (C+ [ ^l^-i;( X )R( X J £ (s)) 2 dx] ^-ds. 
JO ' Js=0 V JR d £ J K 



According to (|20j) . we deduce that 



I I n e (s,x) ^ R ^ x j e ^f dx dg < C + CeT ^ 
Js=0 JR d £ 

and the result (fl~5j) follows. 
The result on g £ (t) = J n £ (t,x)dx follows similar lines. We write 

^6e(t) = \ J n £ (t,x)R(x,I £ (t))dx := K e {t), 



and 

d 



dj KS) = Jn £ AR(x,I £ (t)) + J^R(x,I £ (t)) 2 + J ^R( x ,I £ (t)) ±I e Qk) 

> -K 3 + ^J £ (t) (21) 

>-K 1 -K 3 + ^(j £ (0))_e- K ^ 
(after using the inequality (fl~9j) for (jT^i))). This implies 

K £ (t) > /C e (0) - (K! +K 3 )t- (J £ (0))_. 

The result on —g £ (t) follows. □ 

4 Constrained Hamilton-Jacobi equation, uniqueness 

As already mentioned earlier, in the limit e — > 0, the solution n £ to ([I]) converges weakly to a measure 
n £ L°°(lR + ; L 1 (M (i )) (see see Remark I3.4D . In this section we give a general theory for describing 
properties of n. The statements of Theorem 12.21 follow from the present analysis (and consequences 
in next section. 

We expect that n £ concentrates as Dirac masses 

n £ (t,x) ->> n(t,x) = £;(t)<5(x - £$(*)). 

i 

This weak limit can be described more accurately through the phase function ip £ defined by 

n £ (t,x) = e^ t ^ £ , (22) 

just as the Dirac mass at is well approximated by the gaussian e~~^ /( 2e ). The description of 

V27re 

ip £ , and its limit, gives information on the measure n. 
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Firstly, we obtain the following equation, equivalent to ([!]), 

' d 

—<p e (t,x) = \Vip e \ 2 + R(x,I E (t))+eA(p E , 
< ot (23) 

tp E (t = 0, x) = ife(x) := elnn £ . 

Following |18^ [6], in the limit, we obtain a viscosity solution to the constrained Hamilton- Jacobi 
equation 

J^(t,aO = |Vv9| 2 + i?(x,I(t)), 
< max p(t,sc) =0, Vt > 0, ( 24 ) 
> p(t = 0,x) =y>°(a;). 

Theorem 4.1 issame (E|)-(G|) (v^e * s a sequence of uniformly bounded function in W 1 ' 00 which 
converges uniformly to ip°. Then, after extraction of a subsequence, (tp E ) e converges locally uniformly to 
a Lipschitz continuous viscosity solution ip to (24\ ). In particular, a.e. in t, supp n(t, •) C {<f(t, •) = 0}. 
This solution is unique when assumptions U0\) and 

R(x,I) = b(x) — d(x)Q(I), with Q(I) > increasing, 

or 

R(x,I) = b(x)Q(I) — d(x), with Q(I) > decreasing, 

Because we expect that the set {<p(t, ■) = 0} is made of isolated points, we indeed expect that n is a 
sum of Dirac masses. In one dimension this follows rigorously from Theorem 15.11 below but heuristi- 
cally maximum points of (p(t, ■) are indeed isolated in any dimension. 

The result was proved in [6] in the case R(x, I) = b{x) — d{x)Q{I) and for an integral operator 
instead of a diffusion. For completeness, we recall the proof and the main new ingredients. One of 
the difficulties being that coefficients, because of I{t), are not discontinuous. This can be handled 
following well established arguments, see 

Proof. Existence. In order to prove it, we show how to pass to the limit in f|23|) . It is worth pointing 
out that we can do it in two ways, either by using the notion of viscosity solutions for equations with a 
^-dependence in time or, since I £ converges to an increasing function, by using the notion of viscosity 
solutions for equations with a discontinuous Hamiltonians. We do it in the second way which has the 
advantage to give a more precise result. 

Since the functions f e are equi-bounded and equi-Lipschitz continuous in space, they are also equi- 
Hdlder continuous in time (See [1] or [23] for proofs of this claim). Therefore the only difficulty is to 
pass to the limit (up to a subsequence) in the term R{x,I £ ). 

More precisely, in order to pass to the limit in this term, we have to compute the quantities (up to 
a subsequence) 

limsup I £ '{ s ) and liminf I £ >(s) . 
s->t s->* 
e'-»0 

Indeed, the rr-dependence in R(x, I £ ) does not cause any problem since the functions x \—* R(x, I) are 
Lipschitz continuous in x and, on the other hand, the function R(x, I) is decreasing in / for any x 
which allows to reduce the computation of these limsup and liminf to I £ i and not R(x,I £ '). 
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To do so, we use Theorem 13. II : if we set /x e := lr e i/2 T i + 2Ce) dt, where, here and below, l[ a fe ] 
denotes the indicator function of [a, b], then the // e 's are bounded (positive) measures on (0,T) for 
any T > 0, if e is small enough. Therefore we can extract a subsequence such that (/v)e' converges 
weakly to some measure /x. 

If < to < T is a point such that //({to}) = 0, then we can also assume that I £ '(to) converges to 
some a£l and we set 

I(t) = a + I dfi(s) if t > to an d /(t) = a — dfi(s) if t < to • 

Obviously, the function / is an increasing function on (0,T). 
We are going to show that 

limsup I £ '(s) < I(t+) and liminf I £ >(s) > I{t— ) . 

We do it for t > to, the other case being treated analogously. 
For s close to t, we have 

4'( s ) = " + / ^(t) • 

J(t ,i) 

If C : [0, T] — > M is a continuous function such that £(t) > lr t0it i(r) on [0,T], we have 

I e >(s)<a+ [ C(r)d^(r) , 

J(0,T) 

and passing to the limit in this inequality, we obtain 



limsup I £ >(s) < a + / ((T)dfi(r) , 
7(0,T) 



e'->0 

and we conclude by applying Lebesgue's dominated convergence Theorem to a sequence (Ck)k which 
converges to l[t 0) t] in a suitable way : namely, we choose a sequence of functions such that < < 1 
on [0, T], Cfc = 1 on [to,t] for any k and which converge pointwise to outside the interval [to,t]. 

The proof for the liminf is the same approximating, this time, the indicator function from below. 
It is worth pointing out that the property //({to}) = avoids here (and above) discussions on the 
behavior of the sequence at this point. We find that this liminf is larger than a + f,^ t * ^(r)cZ/i(r), 
which is exactly I{t—). □ 

Proof. Uniqueness. Consider for instance the second case. We consider again the function 

V{t,x) = ip(t,x) - 6(x)£(t), £(t) = f Q(l(s))ds. 

It satisfies 

dt 







-m, x) = -d(x) + |V(* + 6(x)E(t)) | 2 . 
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On the one hand, for two different solutions with the same initial data, we define in this way two 
functions and ^ 2 . Using the viscosity criteria, we have, at the point xq where the maximum is 
achieved 

^||*i-*2(t)||oo = |V(^ 1 (xo)+6(x )S 1 (t))| 2 - |V(^ 2 (x ) + 6(x )S 2 (t))| 2 



= (V*i(x ) + V6(xo)S!(t) + V* 2 (x ) + V6(x )S 2 (t)).(V6(x )S 1 (t) - V6(x )S 2 (t)) 

< C|V^o)Si(t)-V6(x )S 2 (t)| 

< a||V6||oo - S 2 (*)l (25) 
On the other hand, we also have, considering the point Xi where max^g^d (fi(t,x) is achieved 

= max (pi(t,x)— max</? 2 (t, x) 

xeR d x£R d 

< (fi(t, xi) - ip 2 (t, Xl) 

< 6(xi)[Ei(t) - E 2 (t)] + *i(t,xi) - * 2 (t,xi) 

< 6(si)[Ei(t)-S 2 (t)]+niax[*i(t,x)-*2(t,s)]. (26) 



Changing the indices 1 and 2, one of the numbers b(xi)[T,i(t) — S 2 (i)] or &(x 2 )[E 2 (i) — is 
negative and thus we have, with the notation ()10p . 



6 m |(Ei(t) - £ 2 (t)| < max |*i(*,x) - * 2 (i,x)|. (27) 
Together, the two above inequalities yield 

^11(^1 - *2)(*)||oo < C||(^! - ^X^Hoo- (28) 

Uniqueness follows. □ 

The uniqueness result can be obtained under a slightly more general assumption on R, namely: 
for any functions Ii,/ 2 6 L°°(0,T) and t G (0,T), the function x i— > J Q * [i?(x,Ii(s)) — i?(x,/ 2 (s))] ds 



rt 

has a constant sign for all x € R d and there exists a constant K > such that 



max 



o 



[A^MiOO) -D x R(x,I 2 (s))] ds 



< K min 



[R(x,h(s))-R(x,I 2 (s))]ds 



(29) 



The result can be stated as follows 



Theorem 4.2 Assume that cp° e W 1 '°°(R d ) and that tffl-tffy, (E2P and the above sign condition hold. 
Then the problem \2J$ has at most one Lipschitz continuous solution. 



The two conditions appearing in assumption (|29j) are easily checkable in the cases we emphasize 
in Theorem 14.11 but (a priori) they are not satisfied for general R of the form R(x,I) = b(x)Qi(I) — 
d(x)Q2(I), even if ([9]), (TL0j) hold true. Uniqueness for such general R is still an open problem. 

Proof. The proof follows readily the uniqueness proof of Theorem 14.11 We consider two solutions 
(cpi, Ii) and {ip2, 12) and for i = 1, 2, we set 

^>i(t,x) = ifi(t,x) - R(x,Ii(s))ds. 
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The first part of the proof leading to (|25p yields now 



- * 2 (t)||oo < Kmax 
at x pR d 



[D x R(x,h(s)) - D x R(x,I 2 (s))} ds 
and the second part (cf. (]26p -( j2Tj) ) provides us with the inequalities 

0< I [RixxJ^s)) - R(x 1 ,I 2 {s))]ds + max[^ 1 (t,x)-^ 2 (t,x)], 

JO x€R d 

0< / \R(x 2 ,I 2 (s)) -R(x 2 ,h(s))] ds + max [* 2 (t,x) - Vi(t,x)]. 
Jo xm d 



and then 



mm 



[R(x,h(s)) -R{x,I 2 (s))}ds 



< max \^i(t,x) - ^ 2 (t,x) 

xeR d 



(this is where it is important that x i— > [R(x, h(s)) — R(x, I 2 (s))] ds does not change sign). Putting 

together the two inequalities for "^ll^i ~~ ^2(£)||oo and ll^i — ^ / 2(*) ||oo and using assumption (f29l) leads 
to ([28]) and uniqueness follows. □ 



Remark 4.3 In [6], the uniqueness proof is given in full details for merely continuous solutions. 
Lipschitz continuity of ip is used here to simplify the arguments (an write them a.e.) but everything 
can be understood in viscosity sense. 



5 Structure properties 

We now derive consequences of the constrained Hamilton-Jacobi equation in terms the concentration 
points. 

Theorem 5.1 Assume |5]). For any tp° G W l, °°(M. ), the solution to \2J$ is semi-convex in x for 
t > 0, i.e. for any t > 0, there exist a C sc {t) such that, for any £ £ M. d , |£| = 1, we have 

d 2 

-Q^Ve > -C sc (t), i = l,...,d. 

Consequently, at a maximum point x{t) of tp(t, ■) in the variable x, ip(t, ■) is differ entiable in x (but 
maybe not int) and we have 

V(p(t,x(t)) = 0, 
additionally for all Lebesgue points of I(t) we have 

R(x(t),I(t)) = 0. (30) 
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The following proof does not use the BV property of I(t), (only the strong limit I e — > I is used 
to pass to the limit in the Hamilton- Jacobi equation), and thus might be useful for extensions of the 
model where BV estimates are not proved as in section El 

Proof. (Proof of semi-convexity) We justify the semi-convexity directly on the diffusive equation and 

d(p £ d 2 (p £ 
use the notation pt := -— , (pee := , We have 

at, a£ 

d 

—ip^ = 2Vp £ • + R^x, I £ (t)) + eAtpt, 

d 

= 2Vp £ • Vy> K + 2|V^| 2 + %(x, I £ (t)) + eAiptf. (31) 
But \Vip^\ > \<P£z\ because tp^ = Vtpg • £, therefore the function w := (p^ satisfies 

—w > 2V(p £ ■ Vw + 2w 2 - R + eAw 
ot 

with R an upper bound on R^(x, I e (t)). The semi-convexity result (in fact, the semi-convex regular- 
izing effect) follows from suitable approximation arguments (for the initial data) and the comparison 
with the (absolute) subsolution given by the solution of the O.D.E. y = 2y 2 — R, y(0) = — oo. 

(Proof of V(p(t,x(t)) = 0) We prove this equality and a related preliminary result. As we have 
shown above, the function <p is semi-convex and, by classical properties of semi-convex functions, ip is 
differ entiable at maximum points. Therefore 

Vip(t,x(t)) = 0. 

Moreover, it is standard that for any sequence {tk,Xk) of x-differentiability points of (p which converges 
to (t,x(t)), we also have 

Vtp(t k ,x k ) -> 0. 



From this property, we deduce that, for h, r > 0, h, r — ► 0, 

rt+h fx(t)+r 2 ft rS(t)+r 

\V<p(s,y)\ 2 dsdy , — \ 

c(t)-r rh Jt-hJx(t)-r 

Indeed, considering the first limit for instance, a straightforward change of variable yields 

-1 /•! 



— / \V<p(s,y)\ 2 dsdy , — / / \V V (s,y)\ 2 dsdy^0. (32) 

rh J t Jx(t)-r rh Jt-h Jx(t)-r 



J J \Vtp(t + hr,x(t) +re)\ 2 drde , 



and it suffices to apply Lebesgue's dominated convergence Theorem to this integral, using the above 
property together with the (local) Lipschitz continuity of (p. 

(Proof of R(t,x(t)) < 0) We come back to the equation which holds almost everywhere, and we 
first integrate it on rectangles (t, t + h) x (x(t) — r,x(t) + r). We obtain 

rx(t)+r rt+h px(t)+r rt+h rx(t)+r 

/ (<f(t + h,y) - <p(t,y)dsdy / R(y , I (s)) dsdy + / \Vp(s,y)\ 2 dsdy . 

Jx(t)-r Jt Jx(t)-r Jt Jx{t)-r 
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But, by the semi-convexity of <p in x, we have, for y £ (x(t) — r, x(t) + r) 

0><p(t,y) > ip(t,x(t)) - -^\y - x(t)\ 2 = 0(r 2 ) , 



while we have also <p(t + h,y) < 0. Using these two properties in the above equality, we deduce 

i rt+h rx(t)+r i rt+h rx(t)+r i 

I / R(y,I(s))dsdy + — / | V<p(s, y)\ 2 dsdy < — 0(r 2 

t Jx(t)-r rhj t J S (t)-r rh 



rh 

Therefore, we obtain 



i rt+h rx(t)+r i 

/ / R{y,I(s))dsdy< — 0{r 2 ) 

Jt Jx(t)-r rn 

and letting r, h tend to zero with r < /i, we conclude that at any Lebesgue point of / we have 

R(x(t),I{t)) < . 

(Proof of R{t,x{tj) > 0) To obtain the opposite inequality, we integrate on the rectangle (t — h,t) x 
(x(t) — r,x(t) + r). The left-hand side is now 

(<p(t, y) - fit - h, y)dsdy > / (p(t, y)dsdy 

x(t)—r Jx(t)—r 

and, by the afore mentioned argument, it is larger than 0(r 2 ) and we are lead to 

I rt rx(t)+r ^ rt r-x(t)+r ^ Q(r) 



rh 



rl rx(t)+r ^ rt rx(t)+r 

/ / R(y,I(s))dsdy + — / \V<p(s, y)\ 2 dsdy > 

Jt-h Jx(t)-r rn Jt-h Jx(t)-r 



r h Jt-h Jx(t)-r r h Jt-h Jx(t)-r ^ 

Again we let r, h tend to zero with r <S h and using (|32p . we conclude that, at any Lebesgue point of 
/, we have 

Q(t,x(t)) > . 

□ 

Let us recall that because I(t) is BV, it is continuous away from a countable set of discontinuity 
points. Hence, we can give another proof of the relation (|30[) . Let to a continuity point of I{t) and let 
x{to) be a maximum point of ip(to). 

First, we use the viscosity subsolution criteria at such a point (to,x(to)), testing against the test 
function 0. We find 

0<R{x(t ),I(t ))+0. (33) 

To get the other inequality , we integrate in time the equation (I24p on (to, to + h) (h > 0) at the point 
x = x(to) and find 

> ip(to + h, x(t )) > \ I R(x(t ),I(t + s))ds. 
"> n J S=Q 

Because to is a continuity point of /(to), we find 

Q ^ liminf y(to + / t ,x(to))- y (to,x( to )) ^ moim)) {M) 

h^0+ h 

Clearly, ^ and (JM} give the relation ffHTF]). n 
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6 Continuity and discontinuity of I(t) and x(t) 

So far our results have shown that the Lagrange multiplier I{t) is BV and that the concentration 
point x(t) is unique in one dimension. It is natural to ask whether these quantities are continuous. 
This question is particularly relevant because it is possible to derive formally a differential equation 
on x(t), see pi IH] 

j t x{t) = {-D 2 ip{t,x)Y l V x R{x{t)J{t)). 

In this section, we answer negatively to this question (in general) based on analytical and numerical 
examples. We also give a smallness assumption implying continuity. 



6.1 An analytical counter-example 

Our analytical counterexample follows earlier ideas developed for combustion models, see [221 [29] for 
instance. 

We choose in this section t/j = 1, in other words I(t) = g(t). We first consider the simple equation 



together with the initial data 
and the constraint 



ip t = x — g(t) + \ f x \ 2 in (0, +oo) x R, 

ip(x, 0) = (fo(x) in R , 
max (p(t, x) = for any t . 



(35) 

(36) 
(37) 



In a first step, we start by solving the initial value problem for \& = ip — f g(s)ds with the initial data 
tpo(x) = —x 2 . We look for a solution of the form 



Easy computations yield 



a(t)(x - b(t)) 2 + c(t) . 

a'(t) = 4a 2 (t), a(0) = -l, 

-2a(t)b'(t) = 1, 6(0) = , 

[ c'(t) = b(t), c(0) = . 

Therefore, we have the explicit coefficients and solution 

1 , / 



a(t) 



t 2 t 3 



. bit) = - + 1 2 , cit) = — + — , 
1 + At ' w 2 ' w 4 3' 



, * if t 2 y t z t- 1 

V(t,x) = [x r H h — 

v ' ; l + 4t\ 2 / 4 3 



From this, we deduce that the solution to (I35l)-(l371) is 



^x) = - T ^ I (x- t --t 2 



g(t) 



+ t 2 . 



This is consistent with the property g(t) = x(t) = b(t) of Theorem 12.21 
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In a second step, we remark that, in order to solve the problem on ^ = ip — g(s)ds with initial 
data — (x — a) 2 — 5 (a > 0), we have just (i) to translate the solution obtained with a = 0, (ii) to 
subtract 5 — at to take into account the new right-hand side of the equation and the initial data. 
Therefore this solution is 



If t 

x — a v 

1 + At V 2 



t A 



5 + at + — + 

4 3 



Finally, we solve ([55 |1 - (|37|1 with 

(fo(x) = max(— x 2 , — (x — a) 2 — 5) in 



Since the Hamiltonian of (|35p is concave, the maximum of two solutions is a solution : this is a 
consequence of the Barron- Jensen approach ([9l[T0], see also [2]). Thus the solution of this problem is 



1 



max 



1 + 4t V 2 
with g(t) = \ + t 2 , while 



t 



If t 

\ x — a V 

1 + At I 2 



5 + at] for t < t := - , 
/ a 



1 



max 



1 + At 



t 

x — — — t +o — at, 



If t 

— [ x — a V 

1 + At V 2 



with g(t) =a + ± + t 2 . 



for t > t , 



This shows the phenomena of discontinuity of g(t) and of the Dirac concentration point x{t). 




Figure 1: A numerical example of discontinuous solution in example (f38|) . The abscissae is time. (Left) 
The Lagrange multiplier g(t). (Right) The concentration point x(t) (in fact isovalues of the density 
n e {t,x)). 



6.2 A numerical counter-example 

We give a numerical simulation of a similar phenomena which however does not use a discontinuity 
of the gradient to generate the discontinuities of the Lagrange multiplier and the concentration point. 
We use again tp = 1 (J = g), and 

R(x, g) = ( x -x 2 + 3x 4 )(9 - (1 + x) 3 ) - g. (38) 
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Figure 2: Discontinuity in example (|38|) as in Figure [T] The phase <p(t,x) as a function of x at different 

TIMES AROUND THE DISCONTINUITY TIME. We OBSERVE THAT THE DISCONTINUITIES DEPICTED IN FiGURe[T]cORRESPOND 
TO A SMOOTH TRANSITION ON if. 

The computations depicted in Figures [I] and [2] are performed with < x < 1 and (p°(x) = — \x — 0.05|. 
They show how a smooth change on the phase funtion ip(t, x) along time, can produce a discontinuity 
of the Lagrange multiplier I(t) and the concentration point x(t). A local maximum of ip(t,x) in x 
becomes larger as time goes and becomes the global maximum at the discontinuity time. This follows 
the scenario in the analytical counter-example. 

6.3 Concavity conditions for continuity 

In this section, we are interested in conditions that imply the continuity of the I{t) and x(t). Because a 
discontinuity corresponds to a double maximum of ip(t, x), it is natural to look for conditions implying 
the concavity of the solution ip. This is the purpose of the next 

Theorem 6.1 Assume ([5|j and that R(-,I) is strictly concave, uniformly for bounded I. Then, for 
any ip° £ W 1,oa (R d ), ip° uniformly concave, any solution to \2J$ is strictly concave and thus x(t) is 
continuous. Therefore, in Theorem I5.il the equalities hold everywhere and I(t) is also continuous. 

Proof. Let <p a solution of (I24|) and < a < 1. We consider the function 

X(x, y, t) := a<p(t, x) + (1 - a)ip(t, y) - ip(t, ax + (1 - a)y) . 

The aim is to prove that x( x :V:t) < for any x, y G M d such that x y and for any t. To do so, we 
first show that x is a viscosity subsolution of 

2 2 

Xt < 2K\ Xx + Xy\ + — + 7— ^— + Ra{x, y, t) in R d xR d x (0, +oo) , 
a 1 — a 

where K := \ \DipWoo and 

R a (x,y,t) = max {aR(x, J) + (1 — a)R(y, J) — R(ax + (1 — a)y, J)} . 

I ^1 ^ I l-^l |oo 

We are going to argue below as if the function t *—>■ R(x,I(t)^J were continuous : a rigourous proof 
would consists in approximating / by continuous functions, in proving that the corresponding 
satisfy the above inequality (this is also why we put a max in the definition of R a : to point out the 
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uniformity in I) and finally to pass to the limit. We drop these details here for the sake of simplicity 
in the presentation. 

Consider a test function <p(x, y, t) in the viscosity sense and (x, y,t) is a strict maximum point of 
x(x,y,t) — <p(x,y,t), we look for local maximum points of 

lax + (1 — a)y — z\ 2 
a<p(t, x) + (1 - a)(p(t, y) - ip(t, z) 4>{x, y, t) 

near (x,y,t). By classical arguments, since (x,y,t) is a strict maximum point of x(x,y,i) — (p(x,y,t), 
there exists a sequence (x £ , y £ , , z e , t £ ) of maximum points of this new function which converge to 
(x, y, ax + (1 — a)y, t) as e tends to 0. 

If p := — — — -, standard uniqueness arguments (cf. [15]) provide the existence of 

a,b,c (playing the roles of the i-derivatives of (p at respectively (x e ,t £ ), (y £ ,t £ ), (z £ ,t £ )) such that 

a < lp+ D ^e,ye,t £ ) 2+ ^ 

a 

1 — a ' 
c> H 2 + ^(z £ ,I(t £ )) , 

and with aa + (1 — q)6 — c > cpt(x £ ,y £ , t £ ). 
Combining these inequalities, we obtain, after straightforward computations 



(t>t{x £) y £ ,t £ ) < 2p ■ (D x (j)(x £ , y s , t £ ) + D y <p(x £ ,y £ , t £ ))+ 



D x (j)(x £ ,y £ ,t £ 



D y (j)(x £ ,y £ ,t £ 



1 — a 



+ aR(x £ ,I{t £ )) + (1 - a)R(y £ , I(t £ )) - R(z £ , I{t £ )) . 



To conclude, we first remark that \p\ < \ \DipWoo since p is in the x-subdifferential of ip at x, and then 
let e tends to 0. Since z £ — > ax + (1 — a)y, we have 

ai?(x £ , I(t e ))+(l-a)R(y e , I(t e )) -R(z £ ,I(t £ )) -» ai?(x, /(t))+(l-a)i?(y, /(f)) -i2(ax+(l-a)y, /(t)) 
and we are done. 

In order to prove that x( x > 2/j t) < f° r an y x iV ^ ^ d such that x ^ y and t > 0, we first remark 
that x{ x iDit) — f° r an Y x, y G M d as a consequence of a standard comparison result : indeed since 
R a (x,y,t) < (because R(-,I) is concave for any I), x is a subsolution of the P. D. E. 

u? w 2 

w t = 2K\w x + u>„| + — + — v — in R d x R rf x (0, +oo) , 
a 1 — a 

while is a solution, and x(x, y, 0) < in M rf x M d because is strictly concave; therefore x{ x i Vi^) — 
in R d x R d x (0,+oo). 

On the other hand, if x( x i Hit) = at some point (x, y, t), then (x, y, t) is a maximum point of x ~ 
and the viscosity subsolution inequality reads 

< R a (x,y,t) , 

a contradiction with the strict concavity of R(-,I), uniformly in I for / bounded, which implies 
R a (x,y,t)<0. n 



17 



7 The monomorphic case: more precise statements 

We consider in this section the special case when d = 1 and we assume a monomorphic situation 
according to (|12p which leads us to introduce the following assumption 

There exists a function R{x, I) such that x *— > R(x, I) is strictly monotone for any I, I \—* R(x, I) is 
strictly decreasing for any x and 

R(x, I) = (resp. > 0, < 0) if and only if R(x,I) = (resp. > 0, < 0). (39) 



A typical example of such situation was given in the introduction, see ([8]) with assumptions (l9l)~ (jT0j) . 
When R is of the form R{x, I) = b(x)Q 1 (I) - d(x)Q 2 (I). Then ([39]) reduces to 

b(x) 

x i ^ — - is strictly monotone. (40) 

Indeed, one can choose R(x,I) = { ' } = — ^\ } . This is weaker than assuming b' > 

d(x)Qi(I) d[x) Qi(I) 

and d! < for which R satifies (f39j) . 

Theorem 7.1 ^sswme and that tp° £ PF 1,0O (R d ) satisfies 

max(/9°(x) = and limsup ^ (a;) < . (41) 

K |x|^+oo 

For all t > except perhaps for a countable number, there exists a unique x{t) such that 

tp(t, x(t)) = max cp(t, x) = . 

Moreover 1 1— > is monotone and, with the notations x(t + ) = lim s ^x(s), a;(t _ ) = lim s ^a;(s) ; i/ien 

//-R is increasing in x then x{t + ) the largest maximum point of <p(t, ■) and x{t~) is the smallest 
one. 

2. If R is decreasing in x then x(t + ) the smallest maximum point of ip(t, •) and x(t~) is the largest 
one. 

Proof. Using assumption (|4ip and cone of dependence type properties (recall that ip is Lipschitz 
continuous), we see that, for any T > and t £ (0, T), the maximum of the function x \— > <p(t,x) is 
achieved in a fixed compact subset of R. And, by (|30p . if i is a Lebesgue's point of I (meaning here a 
continuity point of the increasing function I) and x is maximum of cp(-,t), we have 

R(x,I(t)) = . 

Therefore, by (|59"]). 

R(x,I(t)) = 0; 

we deduce that x is unique for such t, we denote it by x(t) and, from ([9|) and (|40|) . we see that 1 1— ► 
has the same monotonicity as i? since I(i) is increasing. 

If tfc | t, passing to the limit in <p(tk, x(tk)) = 0, we see that ip(t, x(t + )) = 0. Choosing a sequence 
(t&)fc of continuity points of x(-) or equivalently of /(•), we have also 

R(x{t + ),I(t + )) = . 
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Now assume that R is strictly increasing in x and that there exists a maximum point x > x(t + ) of 
(p(t, •). Then for y close enough to x and s > t, using that I(t + ) = hm s | t I(s), we have 

%,J(s)) > R{x(t+),I(s)) +7?, 

for some 77 > 0, and therefore, using again (|39l) . for s close enough to i 

fl(y,/(*)) >r?, 

for some fj, which would imply that </?(s, x) > 0, a contradiction. Therefore such x does not exist and 
actually, x(t + ) is the largest maximum point of tp(t, •). 

In the same way, x(t~) is the smallest maximum point of (p(t, •) and the case when R is strictly 
decreasing follows along the same lines. 

8 Several environmental unknowns Ie (t) 

Several possible extensions are possible which are more realistic and give more interesting structures. 
In particular the selection rate R may depend on several integrals and competition may be modeled 
by convolution terms ([231 [TBI I26j). Here we study the former extension when the model contains 
several environmental unknowns. This means that we consider the equation 

^-n £ -sAn £ = —R(x,I £ (t)), x e R d , t > 0, 

m e ( 42 ) 

^ n e (t = 0) = n° £ e L 1 (IR d ), n° > 0, 

and now the environment is described by a vector valued unknown I £ = (/J 1 ^, Ie), with 

j(*)(t) = J \p( { )(x)n £ (t,x)dx, i = l,...,N. (43) 

The t/j^'(x) are still given positive functions (see the precise assumptions below). 

We may have in mind the example when mortality and birth terms depend differently on the total 
population, i.e., N = 2 and 

R(x,I) = b(x)Q 1 (l^(t)) - d(x)Q 2 (lW(t)). 

We may also have in mind the example of the chemostat with several nutrients 

l 



R(x,I) = Y J b l {x)Q l {l {i \t)) - d(x). 
i=i 



Our first goal is to derive the L°° estimates. This uses the following assumptions; there are constants 
ipm < ipM, < J m < Im < 00 such that (vector inequalities should be understood componentwise) 

< Vm < il> < $M < 00, G W 2 '°°(R d ). (44) 

min R(x, I) < 0, when min 1^' > Im, (45) 

Ki<N 
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(46) 



-K < 



a 



R(x,I) < -K- 1 < 0, 



R(-,I)\\w 2 '°°(TBL d ) < K. 



(47) 



01(0 



sup 

I m /2<I<2I M 



Then, an easy adaptation of the proof of Theorem 12.11 yields 

Theorem 8.1 With the assumptions [44\)~(47\) an d I m — Ce 2 < I e (0) < Im + Ce 2 (componentwise), 
we have the uniform bound 



with C a positive constant. 

It is unclear if one can expect BV estimates (and monotonicity results in the limit) on the quantities 
Is • Assuming such a bound, one might follow the arguments of Section O and derive again the 
constrained H.-J. equation (|24p where now I(t) is a N dimensional Lagrange multiplier associated 
with the (single) constraint max x£ ^d ip(t, x) = 0. Of course this rises several fundamental questions 
concerning uniqueness of the solution to this constrained H.-J. equation. Some information is clearly 
lost because the quantities I® are related together in that case when n £ concentrates on less than 
TV Dirac masses (for instance in the monomorphic case). Therefore we face here a fundamental non- 
uniqueness situation. 
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